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ABSTRACT 

We present numerical simulations of the passage of gas through a galactic spi- 
ral shock and the subsequent formation of giant molecular clouds (GMCs), and the 
triggering of star formation. In these simulations, we take account of the observed 
inhomogeneity, or dumpiness, of the pre-shock interstellar medium. As might be ex- 
pected, the spiral shock forms dense clouds while dissipating kinetic energy, producing 
regions that are locally gravitationally bound and collapse to form stars. But the effect 
of the dumpiness of gas as it passes through the shock is to generate chaotic inter- 
nal motions in the gas. The kinematics of these motions are found to agree with the 
observed velocity-dispersion/size relation found in star-forming regions. In contrast 
to the standard picture where continuously driven turbulence generates the density 
inhomogeneities in star-forming clouds, we find here that it is the dumpiness of the 
interstellar gas that produces the chaotic motions as it passes through the spiral shock 
and initiates the star formation process. The velocity dispersion can be understood 
as being due to the random mass loading of clumps as they converge in the spiral 
shock. Within these clouds both the timescale for the decay of these motions, and the 
timescale for forming stars, are comparable to the clouds' dynamical lifetimes. In this 
model there is no need for any internal or external continuous driving mechanism for 
the 'turbulence'. In addition, the coupling of the clouds' internal kinematics to their 
externally triggered formation removes the need for the clouds to be self-gravitating. 
Indeed, while clearly some parts of the clouds are self-gravitating and able to form 
stars, most of the molecular material remains gravitationally unbound. This can pro- 
vide a simple explanation for the low efficiency of star formation. 

Key words: stars: formation — galaxies: ISM — ISM: clouds — galaxies: star clusters 
— galaxies: kinematics and dynamics — open clusters and associations: general. 



1 INTRODUCTION 

Our understanding of the star formation process has recently 
undergone a paradigm shift from the earlier quasi-static pic- 
ture where the environment plays no role, to an awareness 
that star formation is an extremely dynamical process where 
the local environment, and interactions, play a dominant 
role in determining the resultant stellar properties (Hart- 
mann et al. 2001; Larson 2003; Mac Low & Klessen 2004; 
Bonnell et al. 2001). Traditionally, giant molecular clouds 
(GMCs), the sites of star formation, were thought to be 
long-lived entities (ages > 10 dynamical times) in order to 
explain the low star formation rate in our Galaxy (Zuck- 
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erman & Evans 1974; Blitz & Shu 1980; Leisawitz 1990). 
Recently, it has become clear that observational evidence 
from young stellar populations implies a relatively fast star 
formation process whereby clouds appear, form stars and 
then disperse on their local dynamical times of a few million 
years (Elmegreen 2000; Hartmann et al. 2001). One of the 
basic properties of these clouds is that they display internal 
chaotic motions which are highly supersonic. These motions 
are usually interpreted as 'turbulence', and observations sug- 
gest that any driving of the observed turbulence should 
come from the largest scales of the GMC (Brunt 2003). A 
succession of numerical simulations (Mac-Low et al., 1998; 
Padoan & Nordlund, 1999; Ostriker, Gammie & Stone, 1999) 
have demonstrated that such supersonic motions, in the ab- 
sence of some continuous driving mechanism, dissipate on 
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one or two dynamical timescales. The realisation that the 
lifetimes of these clouds are comparatively short then sug- 
gests that the formation mechanism for GMCs is a dynami- 
cal one (Roberts 1969; Shu et al. 1972; Pringle et al. 2001). 
One implication of such a mechanism is that gravitational 
clouds need not control their own evolution in the sense 
that, contrary to the usual assumption, they need not be 
self-gravitating entities (although of course some parts of 
then do need to be in order to form stars), nor in virial 
equilibrium. Thus the internal motions which we see may 
be kinematic, rather than dynamic (Pringle et al. 2001). If 
molecular clouds are globally gravitationally unbound, the 
perceived problems of their low star formation efficiencies 
and supposed long lifetimes are erased (Elmegreen 2000). 
Recent simulations of unbound molecular clouds have shown 
how star formation can proceed while the majority of the 
gas escapes due to its excess kinetic energy (Clark & Bon- 
nell 2004; Clark et al. 2005). 

Star formation has long been known to occur primarily 
in the spiral arms of disc galaxies (Baade 1963) . Spiral arms 
are denoted by the presence of young stars, HII regions, dust 
and giant molecular clouds, all signatures of the star forma- 
tion process (van den Bergh 1964; Schweizer 1976; Bash, 
Green & Peters 1977; Elmegreen & Elmegreen 1983; Rum- 
stay & Kaufman 1983; Ferguson et al. 1998). What is still 
unclear is the exact role of the spiral arms in inducing the 
star formation. Is it simply that the higher surface density 
due to the orbit crossing is sufficient to initiate star forma- 
tion, as in a Schmidt law, or do the spiral arms play a more 
active role? Roberts (1969) suggested that the spiral shock 
that occurs as the gas flows through the potential minima 
triggers the star formation process in spiral galaxies. Shock 
dissipation of excess kinetic energy can result in the forma- 
tion of bound structures which then collapse to form stars. 
This of course can provide the dynamical mechanism we 
need to form the observed molecular clouds. The formation 
of giant molecular clouds in spiral arms can occur in one 
of two ways (e.g. Blitz, & Rosolowsky 2004): through the 
agglomeration of smaller molecular cloudlets as they pass 
through the spiral arms (eg., Pringle et al. 2001) or through 
the formation of molecules in dense atomic gas clouds that 
themselves have been compressed in the spiral arms (Bergin 
et al. 2004, Elmegreen 1993). 

In either formation scenario, we need to explain the 
observed properties of GMCs. GMCs are observed to con- 
tain highly supersonic motions and a wealth of structure 
on all length scales (Larson 1981; Blitz & Williams 1999; 
Elmegreen & Scalo 2004). The supersonic motions are found 
such that the velocity dispersion v varies with length scale 
R according to the relation v oc R 1/2 (Larson 1981; My- 
ers 1984; Miesch & Bally 1993; Heyer & Schloerb 1997; 
Heyer & Brunt 2004). A number of authors, on the basis 
on numerical simulations, argue that it is these supersonic 
motions, maintained by internal or external driving mecha- 
nism, which induce the observed density inhomogeneities in 
the gas (Mac Low & Klessen 2004; Elmegreen & Scalo 2004), 
and that it is therefore the supersonic motions which drive 
star formation. Suggested candidates for an internal driving 
mechanism include feedback from low-mass star formation 
(Silk 1985) even though GMCs with and without star for- 
mation have similar kinematic properties (Williams, Blitz 
& McKee 2000) . External candidates include supernova and 



superbubbles (Wada & Norman 2001; Elemegreen & Scalo 
2004). Although there is sufficient energy in these events 
to explain the kinematics of the ISM and they can poten- 
tially generate the correct velocity dispersion sizescale rela- 
tion (Kornreich & Scalo 2000), what is unclear is their rele- 
vance to star formation, and in particular to the triggering 
of star formation in spiral arms. 

There are, as mentioned above, both observational 
and theoretical reasons for supposing that molecular cloud 
lifetimes are quite short (Elmegreen 2000; Pringle, Allen 
& Lubow 2001; Hartmann, Ballesteros-Paredes & Bergin, 
2001), being at most a few dynamical, or kinematic (cross- 
ing) timescales. If this is the case, then a continuous source 
of turbulence is not required. Indeed, all that is required is 
some initial impulsive input of energy. Moreover, in this case, 
the motions do not need to be 'turbulent' in the usual sense, 
as there is barely time for energy to work its way through 
a cascade of eddies. What is the physical mechanism which 
gives rise to such an impulse? Since star formation occurs in 
spiral arms, the obvious place to look for such a mechanism 
is in the spiral arms itself. 

In view of all this, we suggest here an alternative ori- 
gin for the observed chaotic supersonic motions. It is well 
known that the interstellar medium is highly structured on 
all scales (Lauroesch et al. 2000; Dickey & Lockman 1990; 
Elmegreen & Scalo 2004). The result of passing a shock 
through such clumpy gas is to give rise to internal motions. 
The magnitude of the motions so generated depends on the 
degree of inhomogeneity in the pre-shock medium. We pro- 
pose here that the passage of the clumpy interstellar medium 
through a galactic spiral shock not only produces the dense 
environment in which molecular clouds form (Cowie 1981; 
Elmegreen 1979, 1989, 1991), but also gives rise at the same 
time to their supersonic internal motions. We stress again 
that these motions are not really driven turbulence in the 
usual sense of energy being passed from large eddies down 
through a spectrum of smaller ones (Tennekes & Lumley 
1972). Rather, the passage of a shock through the clumpy 
pre-shock medium leads to supersonic internal motions be- 
ing generated at all scales simultaneously. There have been 
several numerical studies of the effects of spiral arms on gas 
flows. Some assume that the pre-shock ISM is smooth and 
concentrate on the generation of instabilities (e.g. Kim & 
Ostriker 2002) . Others consider global effects of the galactic 
flows (Wada & Koda 2004; Gittens & Clarke 2004). Here 
we present numerical simulations capable of following both 
the galactic-scale gas dynamical response to the passage of a 
spiral arm and the small-scale triggering of star formation, 
but taking account of the inhomogeneity of the pre-shock 
ISM 



2 CALCULATIONS 

We use high resolution Smoothed Particle Hydrodynamics 
(SPH) (Monaghan 1992, Benz et al. 1991) simulations to fol- 
low the passage of interstellar gas through a spiral arm. The 
code has variable smoothing lengths in time and in space 
and solves for the self-gravity of the gas, when included, 
using a tree-code (Benz et al. 1991). 
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2.1 Preliminary simulations of clumpy shocks 

In order to demonstrate the basic idea, we first report on 
calculations which compare uniform and clumpy gas clouds 
passing through a shock induced by a one-dimensional si- 
nusoidal potential. These demonstrate the role played by 
the dumpiness of the pre-shock medium in generating the 
post-shock velocity dispersion. These three-dimensional sim- 
ulations follow the evolution of a non-self gravitating gas 
through a one-dimensional potential of the form 

$ = Acos(kx), (1) 

with k — n/2 and A — 80 such that the potential has a min- 
imum at x = 2. Thus, the dimensionless velocity associated 
with the potential, that is the velocity acquired by falling 
from the peak to the trough, is V po t = v2 x A ~ 12.5V, 
where V is the code units for velocity. This potential was 
chosen to mimic the effect of the passage through a spiral 
arm, but in a more controlled manner. The cloud, repre- 
sented by 2 x 10 5 SPH particles, initially has length 3 units, 
width 2 units and height 2 units. It starts off centred on the 
peak of the potential at x = 0. The internal sound speed of 
the gas is 0.3V, and all particles are given an initial velocity 
of 15V in the :r-direction, corresponding to a Mach number 
of 50. The gas is assumed to remain isothermal throughout. 
These velocities are chosen so that the material shocks itself 
as it climbs out of the trough of the potential in a manner 
analgous to gas passing through the gravitational potential 
of a galactic spiral arm (see below). We carry out three sim- 
ulations, one in which the pre-shock density is smooth, and 
two in which it is clumpy. The two clumpy runs differ in 
the size of the clumps, with radii of 0.1 and 0.2 respectively. 
Both clumpy simulations are constructed of 1000 individual 
clumps of 200 particles each. The clumps are bounded by 
an external pressure term. 

In Figure 1 we show the column density distributions of 
the smooth and small clump simulations as they leave the 
potential minimum and start to shock (t = 0.15) and also 
at a later time (t = 0.25) when the shock is fully developed. 
The Mach number of the shock is w 30. We compute the 
post-shock velocity dispersions from the velocity component 
parallel to the pre-shock velocity. We average the velocity 
dispersion over numerous distinct regions of a given size, all 
centred on SPH particles with high density, and considering 
only particles with large densities indicating their location 
inside the shock. We then repeat this process on many dif- 
ferent sizes to determine how the induced velocity dispersion 
depends on observed size-scale. We show these in Figure 2 
for the three shock simulations at a time (t = 0.25) when 
the shock is fully developed, corresponding to the latter im- 
ages in Figure 1. We see that the uniform shock produces a 
constant, sub-sonic velocity dispersion on all length-scales, 
except for those comparable to the individual particle sepa- 
ration. This is as expected, and demonstrates that we have 
enough resolution to overcome the usual SPH y/N particle 
noise. This velocity dispersion is related to the SPH handling 
of shocks which leave a residual velocity dispersion as it is 
a particle based method which is thus inherently clumpy on 
the smallest scales. 

In contrast, and as expected, the clumpy shock gives rise 
to a significant supersonic velocity dispersion in the post- 
shock gas. The post-shock velocity of a small element of gas 




Figure 1. The column density distribution for the linear shock 
tests are shown for the uniform (top) and clumpy (bottom, with 
r = 0.1 clumps) gas distributions. The left panels (1.2 < x < 
3.2, — 1 < y < 1) show the gas as it leaves the potential minimum 
and begins to shock (f = 0.15). The right panels (2.6 < X < 
4.6,-1 < y < 1) show the gas distribution when the shock is 
fully developed (t = 0.25). The logarithmic column density levels 
range from a minimum of 90 to a maximum of 9 X 10 6 in units of 
particles per unit area. 



depends on the amount of mass it encounters, i.e., its mass 
loading. Conservation of momentum and the variation in the 
amount of mass loading due to the clumpy mass distribution 
thus lead to a velocity dispersion in the gas. 

Within a region of size-scale less than structures in the 
gas, each individual parcel of gas encounters a similar col- 
umn density in the shock. The uniform mass-loading then 
results in a zero or very small velocity dispersion. As the size- 
scale of the region increases, the gas samples a variation in 
column densities due to clumps. Then parcels of gas in the 
region will encounter differing amounts of mass-loading as 
they enter the shock and thus differing decelerations. The 
resulting velocity dispersion thus increases with the size of 
a region. 

This velocity dispersion increases with the size of a re- 
gion as larger regions then encounter a larger range of mass- 
loading in the shock. In other words, larger regions include 
more random samplings from the full range of mass loadings 
and thus post-shock velocities, resulting in a velocity disper- 
sion which increases as R 1/2 . The necessary condition being 
that the clumps have a significant probability of collisions in 
the shock. Thus, from the smallest sizes where the velocity 
dispersion is subsonic, the velocity dispersion increases up to 
sizes of twice the clump diameter, and in so doing becomes 
highly supersonic. On larger scales, the velocity dispersion 
is constant as the mass distribution is effectively uniform 
on these scales. The simulation with smaller clumps (larger 
densities) produces a higher velocity dispersion and a larger 
gradient than does the simulation with larger clumps, but it 
saturates at a smaller size scale. The lower velocity disper- 
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Figure 2. The one-dimensional velocity dispersion of the post- 
shock gas as it passes through the linear potential is plotted 
against size-scale for a uniform (long-dashed) and two clumpy 
shocks, with clumps of radii 0.1 (solid) and 0.2 (short dashed). 
The velocity is in units where the sound speed is 0.3 and the 
distance in units where the spacing between minima in the po- 
tential is 4 and the initial configuration has a maximum spa- 
tial extent of 3. The dotted-line illustrates the velocity dispersion 
size-scale («di S p <x i? 1//2 ) relation deduced for molecular clouds 
(Larson 1981; Heyer & Brunt 2004). 



sion in the larger clump simulation is due to the increased 
filling factor that reduces the dispersion in mass loading that 
occurs in the shock. 



2.2 Flow in a galactic potential 

The galactic potential we use is a combination of a two- 
armed spiral potential taken from Cox & Gomex (2002) of 
amplitude n = 1 atom cm" 3 , and a logarithmic potential 
(eg. Binney & Tremaine 1987) that provides a flat rotation 
curve with v = 200 km/s. The spiral arms rotate as a fixed 
potential with pattern speed 2 x 10 -8 rad yr _1 . In addition 
to the main scientific simulations, we perform a number of 
control simulations in which self-gravity, or alternatively the 
fluid shock, are turned off in order to assess the physical im- 
portance of each in our results. The simulations were carried 
out on the United Kingdom's Astrophysical Fluids Facility 
(UKAFF), a 128 CPU SGI Origin 3000 supercomputer. 



2.2.1 Test particle simulation and the SPH initial 
conditions 

The initial positions of particles in the spiral hydrodynamic 
simulations were taken from a test-particle simulation of or- 
bits in the galactic potential. The particles were initially 
placed on tangential orbits with kinetic energies adapted 
from circular orbits so as to be consistent with the spiral 
potential together with an additional 5 % dispersion. They 
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1.5 x 10 5 


2.5 x 10 6 
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Figure 3. The distribution of mass as a function of gas den- 
sity (points) is plotted before passing through the spiral shock 
(crosses) and once star formation has been triggered (filled pen- 
tagons) for Simulation A. The corresponding cumulative distri- 
butions are shown as the solid and dashed lines. 



were then allowed to evolve over several orbits. The loca- 
tion of a significant density peak of size 100 pc in the spiral 
arm, at a galactic radius of « 4 kpc, was chosen to iden- 
tify « 48000 test particles of interest. The earlier positions 
of these particles, before entering the spiral arm, were then 
used to establish the initial conditions of the SPH simula- 
tions. These test particles were given an additional 5 km/s 
Gaussian vertical velocity dispersion and then each subdi- 
vided into 9 SPH particles to provide the required resolution 
of 4.3 x 10 s SPH particles. These initial conditions were used 
for our simulations. There are three scientific simulations 
which we discuss here, listed in Table 1. The first, Simu- 
lation A, uses the full initial conditions but with particle 
masses chosen to provide a mean surface density of E = 1.0 
Mq pc -2 and a total mass of 10 6 M Q . The second, cumu- 
lation B, is identical except for a lower surface density of 
E = 0.1 M Q pc" 2 and total mass of 10 5 M Q . The third, 
Simulation C, uses only the central region of the above ini- 
tial conditions and then each particle is subdivided further, 
providing a higher resolution simulation of 2.5 x 10 6 SPH 
particles and a mass of 1.55 x 10 5 M and E = 1.0 Mq 
pc -2 . All subdivided particles had the same kinematics as 
the original particles. 

The gas is initially clumpy on scales of several pc with 
a mean surface density of IMq pc -2 , a volume-averaged 
density of 2.5 x 10 -3 Mqpc -3 , a median particle density of 
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10 _2 Mqpc -3 and peak densities in the clumps of w O.IMq 
pc~ 3 (See Figure 3). This dumpiness corresponds roughly to 
the properties of the interarm gas (Dickey & Lockman 1990). 
Self-gravity is initially unimportant due to the large shear 
velocities. For simplicity, the gas temperature was taken to 
be 100 K and remains isothermal throughout the simula- 
tion. This corresponds to a sound speed in molecular gas of 
0.6 km/s. We consider that the gas is either in pre-existing 
molecular cloudlets (Pringle et al. 2001) or that molecular 
formation is sufficiently rapid as to occur in the shock com- 
pression (Bergin et al. 2004). 

2.2.2 Sink particles and star formation 

Star formation is modeled by the inclusion of sink-particles 
(Bate, Bonnell & Price 1995) that interact only through self- 
gravity and through accretion of any gas particles that fall 
within their sink radii. Sink-particle creation occurs when 
dense clumps (p > 10~ 18 gcm~ 3 , 1.4 x 10 3 M Q pc -3 ) of self- 
gravitating gas of size < 0.5 parsecs are collapsing (sub- 
virial). It should be noted that while the internal velocity 
dispersion of these regions is probably underestimated as 
the particles all lie inside one to two SPH kernels, the cho- 
sen temperature of 100K mimics the internal support ex- 
pected for lower temperature gas with supersonic motions. 
This ensures that the star formation modeled here is real. 
The sink-particles have initial (minimum) masses of order 25 
to 100 Mq and therefore cannot be considered as individual 
stars but can be thought of instead as stellar clusters. Ac- 
cretion quickly increases the mass of these 'clusters' to 10 2 
to 10 3 Mq) with maximum masses of m w 10 4 Mq. The 
detailed computation of star formation within these clusters 
is beyond the scope of these investigations. 



3 TRIGGERING OF STAR FORMATION IN 
THE SPIRAL SHOCK 

The evolution of the gas over 34 million years as it passes 
through the spiral potential is shown in Figure 4 for Simu- 
lation A with an initial surface density of E = IMq pc~ 2 . 
The initially clumpy, low density gas (p w 0.01 M Q pc~ 3 , 
10~ 24 g cm~ 3 ) is compressed by the spiral shock as it leaves 
the minimum of the potential. The shock forms some very 
dense regions, which due to the accompanying dissipation of 
kinetic energy and the increased importance of self-gravity, 
can become gravitationally bound and thus collapse to form 
regions of star formation. These regions, in excess of 10 3 
Mope -3 are replaced by sink particles. Further accretion 
onto the sink particles raises their masses to that of typi- 
cal stellar clusters (10 2 to 10 4 Mq). Star formation occurs 
within 2 x 10 6 years after sufficient densities are reached for 
the gas to be recognised as a molecular cloud (10~ 22 to 10~ 21 
g cm -3 ). The total spiral arm passage lasts for ~ 2 x 10 7 
years. The gas remains globally unbound throughout the 
simulation and re-expands in the post-shock region. 

In our simulations, the total fraction of gas turned into 
stars varies from between 5 and 30 % depending on the sim- 
ulation, total gas mass and surface density. In reality this 
is likely to be an upper limit since there are other mecha- 
nisms not included in our computations such as the effects 
of winds and radiation from massive stars and supernova 



explosions. We find that the total lifetime of the molecular 
clouds, measured from the time that significant mass attains 
typical GMC densities of 10~ 22 to 10 -21 g cm -3 , to be of 
order 10 7 years, or a few dynamical times {td yn ~ 4 x 10 6 
years) . 

The shock dissipation in the spiral arm is essential to 
trigger the star formation process. We have also run simu- 
lations that exclude the shock (with no SPH artificial vis- 
cosity) but include self-gravity; in these we find that there 
is no induced star formation even with orbit crowding and 
the increased importance of self-gravity in the spiral arms. 
Similarly, in Simulation B which is a full fluid simulation 
with but with a lower initial surface density of E = O.IMq 
pc~ , only limited star formation occurs, producing two 'star 
clusters' before the gas leaves the spiral arm and disperses. 
Triggering of star formation thus requires both the dissipa- 
tion of the excess kinetic energy in the shearing flow and a 
critical threshold in surface density in order for self-gravity 
to become important. 

In addition to forming dense clouds in which star for- 
mation occurs, the spiral shock forms structures which re- 
semble the observed structures in GMCs. Detailed images 
of the GMCs formed in a spiral shock are shown in Figure 5 
from Simulation C, with the gas density as in Simulation 
A, but at higher resolution simulation with 1.5 x 1O 5 M0 of 
gas, E = 1M Q pc~ 2 , modelled with 2.5 x 10 6 SPH particles. 
The projections of the dense gas along a length of 50 pc, 
viewed from within and from above the disc of the galaxy, 
are taken just after star formation has been initiated. We see 
significant amounts of substructure as the dense gas breaks 
up into many components of several parsecs in size and sep- 
aration. The dense regions have column densities typical to 
GMCs. This structure can be understood as being due to 
the combination of the clumpy nature of the pre-shock gas 
along with Rayleigh- Taylor like instabilities which occur as 
gas continues to flow into the shock. 



4 THE GENERATION OF THE INTERNAL 
VELOCITY DISPERSION 

One of the most significant properties of GMCs is their 
supersonic internal motions generally characterised by the 
internal velocity dispersion (Larson 1981; Heyer & Brunt 
2004). Thus, in addition to triggering star formation, any 
mechanism to explain the formation of GMCs must be able 
to explain the origin of the internal velocity dispersion and 
how it depends on the size of the cloud, or region consid- 
ered (Larson 1981; Heyer & Brunt 2004). Here we discuss 
how spiral triggering of star formation produces these kine- 
matics when the pre-shock gas is clumpy. The basic idea is 
that when structure exists in the pre-shocked gas, the stop- 
ping point of a particular clump depends on the density of 
gas with which it interacts, thus some regions will penetrate 
further into the shock, broadening it and leaving it with a 
remnant velocity dispersion in the shock direction. The dif- 
ferent elements of the gas pass through shocks of different 
strengths and different orientations. The amount of mass- 
loading, and hence deceleration that occurs in the shock is 
variable due to the clumps in the pre-shock gas. Smaller 
scale regions in the shock are likely to have more correlated 
momentum injection as well as mass loading and thus small 
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Figure 4. The evolution of cold interstellar gas (m = 10° Mq) through a spiral arm is shown relative to the spiral potential of the galaxy 
(upper left-panel) for Simulation A.. The minimum of the spiral potential is shown as black and the overall galactic potential is not 
shown for clarity. The column densities of the gas vary from 5 X 10 -5 to 0.5 g cm -2 . The gas is shown at times t = 9 X 10 5 , t = 9 X 10 6 , 
t = 1.8 X 10 7 , t = 2.3 X 10 7 , t = 2.8 X 10 7 and t = 3.4 X 10 7 years from the start of the simulation. The 5 additional panels, arranged 
clockwise, show close-ups of the gas at the latter five times. The gas is compressed in the shock and subregions become self-gravitating 
and collapse to form groups of stars (from t = 2.3 X 10 7 years). The cloud produces stars inefficiently as the gas is not globally bound. 



velocity dispersions. Larger regions will have less correlation 
in both the momentum injection and mass loading such that 
there will be a larger dispersion in the post-shock velocity. 
Furthermore, the spiral potential ensures that the gas follows 
converging streams into the spiral shock such that there is a 
significant velocity component present parallel to the shock 
front. Together, this imparts a significant velocity dispersion 
in the post-shocked gas. 

The evolution of the velocity dispersion in Simulation 
A as a function of the size of the region considered is shown 



in Figure 6. The velocity dispersion is calculated from the 
velocity components in the plane of the galaxy within a 
spherical (3-D) region, and then normalised to produce a 
one-dimensional velocity dispersion, as would be observed. 
The velocity dispersion in then averaged over multiple re- 
gions of the same size all of which are centred on dense gas 
> 1(T 21 g cm -3 ). We only use particles whose densities 
exceed a minimum value of (n > 10 cm -3 ; p > 1.7 x 10 -23 
g cm~ 3 ) to ensure we only consider the gas that has entered 
the shock. Lower density gas, even if molecular, is unlikely 
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Figure 5. A 50 parsec region of the dense gas in the spiralshock is shown in two projections, from within the disc of the galaxy (upper 
panel ) and from above the disc (lower panel), taken from a simulation of 3 X 10 s Mq modelled with 2.5 X 10 6 SPH particles. The gas 
has fragmented into many quasi-periodically-spaced dense regions with on-going star formation. The fragmentation occurs due to the 
dumpiness of the pre-shock gas and due to Rayleigh- Taylor like fluid instabilities in the shock. Column densities of the gas vary between 
5 X 10 — 4 and 1 g cm -2 . 



to excite any line emission. The initially low velocity disper- 
sion, of order the sound speed v 3 w 0.6 km/s, increases as 
the gas passes through the spiral shock. At the same time, 
the velocity dispersion increases more on larger scales than 
on smaller scales thus establishing a velocity dispersion size- 
scale relationship which is approximately 

/ R \ 1/2 

^dis P « 1 ( — j km/s. (2) 

By the time that star formation has been triggered, the 
gas contains the characteristic kinematics of giant molecu- 
lar clouds. The passage of the clumpy gas through the shock 
produces both the power-law slope and the magnitude of the 
velocity dispersion. 

In the spiral shock, there is no distinction between 
the velocity dispersion inside a given cloud and the inter- 
cloud velocity dispersion. Thus, presumably the intercloud 
velocity dispersion in the spiral shock will follow the same 
sizescale relation. Once the clouds have passed through the 
shock they expand, with a simultaneous decrease in the 
internal velocity dispersion due to the dissipation of ki- 
netic energy in the internal shocks. The velocity disper- 
sion maintains a K 1 ^ 2 scaling but the magnitude decreases. 
Such clouds, if observed, should somewhat decrease the slope 
of the intercloud velocity dispersion sizescale relation. The 
caveat on this being that such clouds are likely to be increas- 
ing difficult to detect in this epoch due to their decreasing 
surface densities. 

This driving of the internal velocity dispersion due to 
the spiral shock occurs even for simulations where self- 



gravity is not included. We can thus exclude that the ve- 
locity dispersion in our simulations is simply a reflection of 
the virialised nature of self-gravitating clouds. These clouds 
can be far from virialised, in fact completely unbound, and 
still display the same velocity dispersion. Thus we can fully 
attribute the generated velocity dispersion, and its depen- 
dency on the length scale considered, to the passage of the 
clumpy gas through the spiral shock. 



5 DISCUSSION 

The scenario presented here, of a dynamical triggering of 
star formation by the passage of inhomogeneous interstel- 
lar gas through a spiral shock, is consistent with the emer- 
gent viewpoint that star formation is a relatively fast pro- 
cess whereby clouds form out of the ISM, produce stars and 
then disperse, all on their local dynamical, or crossing, times 
(Elmegreen 2000; Hartmann et al. 2001). Few GMCs are ob- 
served without ongoing star formation suggesting that star 
formation must occur soon after the clouds form. Estimates 
based on the ages and distribution of young stellar popu- 
lations imply that star formation itself occurs on a dynam- 
ical timescale (Elmegreen 2000) while the lack of molecu- 
lar gas around older systems suggest that the clouds dis- 
perse after at most a few dynamical timescales (Leisawitz et 
al. 1989). From our simulations we estimate lifetimes of the 
dense clouds of order 10 7 years, or a few dynamical times 
(tdyn ~ 4 x 10 6 years), with pre-star formation lifetimes of 
order 2 x 10 6 years. We find star formation efficiencies in the 
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Figure 6. The velocity dispersion is plotted as a function of 
size at 5 different times during the passage of the gas through a 
spiral shock (Simulation A). The velocity dispersion is plotted at 
4.2 X 10 6 years (crosses), 1.4 X 10 7 years (open triangles), 1.8 X 
10 7 years (open octagans), 2.3 X 10 7 years (filled circles) and 
2.7 X 10 7 years (filled triangles) after the start of the simulation. 
Star formation is initiated at sa 2.3 X 10 7 years. The velocity 
dispersion is calculated from the velocities in the plane of the 
disc, which are then averaged over many different regions of a 
given size that each have central gas densities in excess of 10~ 21 g 
cm~ 3 . The dashed line indicates the Larson relation for molecular 
clouds where a oc R- 1 / 2 (Larson 1981; Heyer & Brunt 2004). 

range 5-30 percent, but note that these are likely to be re- 
duced by including the effects of feedback from young stars 
and/or magnetic fields. 

Our simulations also demonstrate that the observed in- 
ternal kinematics of GMCs can be caused simply by the 
passage of an initially clumpy ISM through the spiral shock. 
There is then no need for any internal or external driving 
mechanism for the 'turbulence' in terms of magnetic insta- 
bilities, gravitational instabilities, stellar outflows, winds or 
supernovae (Mac Low & Klessen, 2004). The velocity disper- 
sion is generated at all scales simultaneously (Brunt 2003) 
and thus no turbulent cascade of energy from larger to 
smaller lengthscales is required. Once generated, these inter- 
nal motions are likely to evolve as in simulations of decaying 
supersonic turbulence (Mac Low et al.1998). A caveat to this 
notion of decaying turbulence is that only 25 per cent of the 
gas in the dense clouds has centre of mass motions with radii 
of curvature of 10 pc or less. Thus the majority of the in- 
ternal motions in the dense clouds do not show evidence for 
significant eddy motions on sizescales of molecular clouds. 

Furthermore, we have seen that GMCs do not need be 
globally gravitationally bound in order to explain their prop- 
erties or their evolution. This relaxation of a generally held 
assumption (that the GMCs are in virial equilibrium) per- 
mits a straightforward explanation for the low star forma- 
tion efficiency in GMCs. The vast majority of the gas is never 
gravitationally bound and thus does not undergo star for- 
mation. Local regions in unbound clouds can dissipate suf- 



ficient kinetic energy to become gravitationally bound and 
thus collapse to form stars (Clark & Bonnell 2004). These 
regions then provide the initial conditions for detailed stud- 
ies of the star formation process (Bate, Bonnell & Bromm 
2003; Bonnell, Bate & Vine 2003). A recent simulation of 
an unbound CMC has shown that star formation efficien- 
cies on the order of 10 per cent are a natural outcome of the 
unbound nature of the clouds (Clark et al. 2005). 

Although our models are explicitly made with two- 
armed grand-design spiral features, the spiral shock origin 
for star formation could equally apply to flocculent spiral 
galaxies. In such systems any gravitational instabilities will 
be sheared into local spiral arms (Elmegreen, Elmegreen & 
Leitner 2003) as occurs in any weakly self-gravitating disc 
(Rice et al.2003). Shocks due to the passage of gas through 
these spiral arms can then trigger star formation much in 
the way described in this paper. 

5.1 Comparison with observations 

Detailed comparisons with observations of spiral galaxies 
(Tilanus & Allen 1993; Allen et al. 1997; Smith et al. 2000; 
Schinnerer et al. 2003) will require global disc calculations 
(Dobbs et al., in preparation) and should ideally include 
more physical processes such as a realistic equation of state, 
feedback from young stars, and the effects of magnetic fields. 
Notwithstanding, there are several predictions from this sce- 
nario which can be tested observationally. 

First, the mechanisms detailed here is of a clumpy 
gas that forms GMCs due to the passage through a spi- 
ral shock. The implication is that GMCs should reside pri- 
marily in spiral arms of galaxies. This is consistent with 
surveys of the outer Galaxy (Heyer et al. 1998) which finds 
an arm-interarm density ratio of GMCs in the Perseus arm 
to be 28:1 (Heyer & Terebey 1998). The overall distribu- 
tion of GMCs in the outer galaxy also reinforces the impres- 
sion that they are confined to spiral arms (Heyer, personal 
communication) 1 . The global distribution of CO gas in the 
Galaxy is also best explained as being predominantly in spi- 
ral arms (Bissantz, Englmaier & Gerhard 2003). 

Second, the internal velocity dispersion is driven pri- 
marily at the shock itself. Therefore the gas should have rel- 
atively low internal velocity dispersions prior to the shock 
that increase rapidly at the shock front and decay monoton- 
ically thereafter. 

Third, the strength and scaling relations of the velocity 
dispersion should be independent of the proximity of young 
stars or other potential sources of turbulence. This appears 
to be the case in the Carina giant molecular cloud (Zhang 
et al. 2001). 

Fourth, the shock-driven origin of the GMCs implies 
that they should have sharper edges on the upwind side 
where fresh material is flowing into the shock, while the 
downwind edge to the GMCs is smoother due to the gas 
dispersal. This can be seen in Figure 5. 

Fifth, a potential implication of this work is that at 
least some component of the the velocity dispersion seen 

1 The presence of GMCs in spiral arms is, at least in our Galaxy, 
not surprising as the spiral arms are generally defined by the 
presence of star formation and hence of star forming GMCs. 
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in the inter-arm gas (Dickey & Lockman 1990) is due to 
this pumping process in the spiral arms. The structure and 
kinematics that result from the spiral shock should then feed 
directly into the next spiral passage. 



6 CONCLUSIONS 

The triggering of star formation by the passage of clumpy 
gas through a spiral arm can explain many of the observed 
properties of star forming regions. The clumpy shock re- 
produces the observed kinematics of GMCs, the so-called 
'Larson' relation. There is no need for any internal driving 
of the quasi-turbulent random motions. The shock forms 
dense structures in the gas which become locally bound and 
collapse to form stars. The clouds are globally unbound and 
thus disperse on timescales of 10 7 years, resulting in rela- 
tively low star formation efficiencies. 
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